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^ \ Abstract 

Compact U(l) theory in 4 dimensions is used to compare the modified iterative and the 
Laplacian fixing to lattice Landau gauge in a controlled setting, since in the Coulomb 
phase the lattice theory must reproduce the perturbative prediction. It turns out that 
on either side of the phase transition clear differences show up and in the Coulomb phase 
. the ability to remove double Dirac sheets proves vital on a small lattice. 

> : 

(N . 

1 Introduction 

One of the main practical advantages one exploits in Monte Carlo studies of lattice QCD is the 
fact that phenomenologically relevant low-energy observables (i.e. masses and matrix elements 
of operators between physical states) may be studied without a need to choose a specific gauge; 
i.e. they may be determined from gauge-invariant correlators, and hence the functional integral 
p 1 may perform a random sampling over all gauge orbits. 

Nonetheless, there are good reasons to address the question how one may transform a 
given ensemble of gauge configurations to a specific gauge; mainly because on such a gauge- 
fixed ensemble gauge- dependent objects like quark- or gluon-propagators may be determined 
too, and the results may be compared with what one gets in perturbation theory or in other 
non-perturbative frameworks, like the Dyson-Schwinger approach. 

An attempt to fix a given configuration to a specific gauge raises several issues. The first 
question to ask is what a lattice analogue of a continuum gauge like the Lorentz condition 

<W = o , (i) 

might be, and the typical answer is to construct a functional on the space of lattice configu- 
rations such that its extremum simplifies in the continuum limit to ([I]), but the choice is, of 
course, not unique. The simplest option is the standard lattice Landau gauge functional 

F \{U}] = ^ E (l + ^ReTr(^(:r))) (2) 

where U^(x+/jl) = Q(xY U^(x) Q^x) is the gauge transformed link variable and V — J2x^i 
but it may be worthwhile to maximize an "improved" functional instead (see e.g. |TJ). The 
second (more technical) question is which procedure shall be used to search for the maximum^] 

1 If there is a unique extremum at all, i.e. if there are no genuine Gribov copies ||. 
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of a given functional like (0), since the latter represents a very high dimensional extremization 
problem, and there is, unfortunately, no known algorithmic solution. 

In view of this problem, Laplacian Landau gauge fixing has been proposed @ - a method 
which does not directly attempt to extremize (0) but rewrites the problem such that, after 
a set of local constraints is traded for a global one, the latter task may be solved exactly. 
The advantage is that the procedure is unique, i.e. the final configuration does not depend 
on the initial gauge. It has led to beautiful results in quenched SU(2), SU(3) and S'[/(4) 
studies [|J, but since there is no firm analytic prediction to compare with, the situation is 
not conclusive. As a result, a good fraction of the research which tries to link the QCD 
confinement-deconfinement phenomenon to the behaviour of specific one- or two-dimensional 
defect structures (e.g. monopoles or vortices) deals with algorithmic issues and there is, to 
date, no agreement on the physics lesson learned. 

In this situation, we decided to address the question in a simpler setting, namely the pure 
compact U(l) gauge theory in 4 dimensions with the standard Wilson action 

S[{U}} = (3 (l-Re(^(x)^(x+/i)^(x+z>) t ^(x) t )) =(3J2 (l-cos(0 D )) (3) 

where (3 = l/g 2 , U^x) = exp(i8 x ^) G U(l) and 9 a = X ^ U = 6 x ^ + 9 x+ ^ v -6 x+i ,^-6 x ^ is the 
plaquette angle. This theory is known to have two phases 0: For /3>/3 c it is in the ordered 
("Coulomb") phase, for (3 < [3 C it is in the disordered ("confined") phase. Contrary to what 
early and some more recent studies indicate ||, the transition seems to be 1st order |7j with 
&~1.01. 

What is important in our context is that the U(l) theory, if quantized on the torus T 4 , 
shows the Gribov ambiguity mentioned above ||. Several numerical studies have tried to pin 
down the effect gauge fixing artefacts have on the mass of the "photon" in either phase. First, 
Coddington, Hey, Mandula and Ogilvie studied the transverse zero-momentum correlator 

C(t) = Y,j:A(^0)A l (y,t) (4) 

x,y i=l 

where they defined the photon field (and we follow them in this choice) through 

A tl (x)=ImU fl (x) = sm(9 Xtfl ) , (5) 

and they found that - at zero momentum - the "photon" behaves as a massless particle for 
(3> (3 C but as a massive one at (3<f3 c 0. Thereby they confirmed an earlier result, gotten in a 
gauge- invariant approach JlO| . Nakamura and Plewnia extended this work to arbitrary spatial 



momentum, i.e. they measured the correlator 

C^(p, t) = £ A M (x, 0)A v (y, t) e- 1 *-^) (6) 
x,y 

for a few low-lying p and observed that the transverse part at /3 = 1. 1 deviates, for large time 
separations, quite drastically from the perturbative prediction, if (^|) is maximized through 
local iterative steps only ||11|| . They also noticed that the result improves, if the local iterative 



procedure is extended to opt for the best one out of several attempts, and this observation 
triggered a long stream of subsequent activities (see below). 

The goal of the present study is to give a fair comparison between Laplacian and iterative 
gauge-fixing in the pure U(l) theory, where both phases may be considered, good statistics is 
easy to obtain and correlators involving all momenta may be determined through fast Fourier 
transform. After a quick review of the gauge-fixing procedures, the comparison is given for the 
Coulomb and confined phases separately, and we end with a short discussion of the relevance 
of our results to the non-Abelian case. 
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2 Iterative versus Laplacian gauge fixing 



We start with a brief outline of the two gauge fixing algorithms which we are going to compare 
w.r.t. the photon correlators implied. 

2.1 Standard and modified iterative Landau gauge fixing 

Modified iterative Landau gauge fixing is an improvement on standard iterative Landau gauge 
fixing. The latter is tantamount to iteratively passing through the lattice and trying to locally 
maximize the gauge functional (fj) which in the case of the compact U(l) theory reads 

F l{ }] = Tr E F M , F x (9) = 1 E ( 2 + cos(^) + cos(^_ AiM )) . (7) 

In this symmetrized version it is obvious that for a given site x the local gauge transformation 

. Ox,[i * 9 x ^-\- c x , $x— p,,fj, ^ 9x—{l,ii Cx (8) 

which maximizes F x (8) may be found analytically. Moreover, doubling the respective c leaves 
the local contribution F x (8) unchanged, thus offering the possibility of an overrelaxation step: 
in practice one sweeps through the lattice in a mixed overrelaxation/maximization strategy 
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In general, this standard procedure gets stuck in a local extremum. If the common belief 
that the "true" gauge copy corresponds to the global maximum |Tj| is correct, an extra effort 



is needed to improve on this situation. Specifically for the case of the compact U(l) theory 
two types of objects are known which prevent the local iterative procedure from reaching the 



true extremum. Bogolubsky et al. have shown |14]] that in the quest for the global maximum 



one needs to suppress both the "double Dirac sheets" (DDS) and the "zero momentum modes" 
(ZMM) of the gauge field, two concepts which we shall explain in detail. 

DDS: Each plaquette angle x>fll , is decomposed into the gauge- invariant electromagnetic 
flux 9 X)IU/ e] — n,7i[ and the discrete gauge- dependent part 2'Kn X:IXV: n Xjlxv = 0, ±1, ±2 [f[5[j . A 
value n Xjlxv = ±1 represents a Dirac string passing through the plaquette and a plaquette with 
n x ,fiv 7^ is called a "Dirac plaquette". A gauge-invariant monopole inside an elementary 
cube is detected via its non-zero flux 2nJ26 n ^u- Since monopole worldlines form closed loops 
on the dual lattice, the Dirac strings normally span cylindrical 2d sheets dual to the Dirac 
plaquettes. However, on a lattice with periodic boundary conditions, it is also possible to form 
a boundary-free Dirac sheet spanning a whole 2-torus. It represents the exceptional case of 
two monopole loops merging on top of each other. A pair of Dirac sheets with opposite flux 
orientation carries zero action, and can be removed by a periodic gauge-transformation JTT|| . It 



is called a "double Dirac sheet" (DDS). In the confining phase, monopole worldlines percolate 



T7j , and so do Dirac strings and sheets. In the Coulomb phase, monopole loops are small, 
Dirac plaquettes are less abundant and DDS are rare. When they appear, they tend to lie 
parallel to the main axes, presumably for entropic reasons, since this orientation requires fewer 
Dirac plaquettes. In any case, a DDS in the ap plane may be detected by counting the total 
number of \w Dirac plaquettes (e puap ^ 0), which must satisfy 

N$> > 2L a L p . (9) 
ZMM: With periodic boundary conditions the "zero momentum modes" of the gauge field 



do not contribute to the action ((^) - the angle of the Polyakov loop is arbitrary, indicating a 
global U(l) symmetry of the action. For gauge configurations representing small fluctuations 
around constant modes it is easy to see that the global extremum of the functional (0) requires 
9^ = 0. The latter condition may be satisfied through a non-periodic transformation 

9 x ,p, -> C Xill = dx,n ~ Cy, mod 27T , c M G ] -7T, 7r[ (11) 

with the specific choice c M = 9^, but in general this is not a gauge-transformation (the Polyakov 
loop is changed). If one insists that (|TT|) shall be a gauge-transformation, the c M need to be 
chosen from the (finite) set 

c,G{0,±—,±— ,...}. (12) 

J-ip ^p 

Had we chosen the other natural definition of the gauge field, A^x) = — ilog(L/^ (x)) = 9 Xjfl 
rather than (^), the transformation ( |i~T]) with c M = 9^ would amount to a shift in A^(x) such 
that its constant Fourier mode is zero. In other words, abandoning the restriction ( |12|) amounts 
to eliminating a constant background field. While this is frequently done (see e.g. [JH]), we 



feel that sticking with the combination (O, |TJ), where the background is suppressed but not 



eliminated, is somewhat closer in spirit to the original extremization task - choosing, out of 



(|I~2j), a non-trivial c M , amounts to switching to another point on the gauge orbit which, in view 
of the global extremum, looks more promising. In practice, one does not expect any sizable 
difference between the two options, since on a large enough lattice the discreteness of c M turns 
into a rather marginal restriction. 

The "modified iterative Landau gauge" (modified ILG) differs from the "standard" ILG 
discussed above in that it tries to eliminate both the DDS and the ZMM as far as possible. 
To implement it, we proceed in bunches of 100 standard sweeps through the lattice followed 
by a ZMM operation ([LI], |T2|) and a check on the DDS condition (|9]) which we sharpen to 

> 1.5L a L p (13) 
to account for additional monopole-antimonopole pairs which may punch holes into the sheet. 



In case ([13]) is met, the whole gauge-fixing procedure is started over from another random 
point on the gauge orbit. This is repeated until a configuration without DDS is obtained. 
Further bunches of 100 extremization/overrelaxation sweeps are applied until the last round 
has increased the functional ([?]) by less than 10 -8 . To prevent the program from getting 
trapped in an infinite loop, a cut-off on the outer loop is introduced which is chosen O(30) 
in the Coulomb phase and O(10) in the confined phase (these figures are tuned to our (3 
values, i.e. they represent an ad-hoc element in the procedure). In such a case the result of the 
random gauge copy which has led to the highest functional value (^) is post-iterated with local 
maximization/overrelaxation steps to make sure that at least a local extremum is reached. 

For reasons which will become clear later on, we decided to include the ZMM suppression 
( PH , |i~2| ) in our "standard" ILG ensemble, too. 

2.2 Laplacian Landau gauge fixing 

The Laplacian Landau gauge fixing method takes its origin from the observation that 
maximizing (|7|) (we specify the discussion to the U(l) case) is equivalent to minimizing 

G[{9}] = Re( Y, 2 - e ie ^ - e~ w '-^) (14) 

x,p 
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and the latter means, more explicitly, to search for a gauge transformation (H) such that 

G[{6}] = Re( e iCx 2e~ iCx - e^e^e"^ - e ^ e -^^ e -^^) (15) 

x,n 

is minimal. Setting <fr x = p x e~ 1Cx , the bracket suggestively looks like the quadratic form 

£4(-A)^ (16) 

x,y 

and the problem (p~5|) would, in fact, be equivalent to searching for the lowest eigenmode 
of the sign flipped covariant Laplacian, were it not for the set of local constraints \<f) x \ = 1 
(or p x = 1) for all x. The Laplacian Landau gauge [[| ignores the latter and constructs the 
transformation (||) from the eigenmode which corresponds to the lowest eigenvalue of (|16|) , i.e. 
it sets c x = — lmlog(0 a ) and keeps only the global constraint J2x I0z| 2 = 1 (because the mode 
may be normalized). 

The Laplacian Landau gauge (LLG) is unique, i.e. the final result does not depend on the 
initial gauge 0, and it is unambiguous in all cases where (a) the lowest eigenmode of the 
covariant Laplacian is non-degenerate and (b) the result satisfies <fr x 7^ for all sites x. For a 
continuous (i.e. interpolating) gauge field A^x), the local ambiguities 0(x)=O are associated 
with Dirac sheets. For details with respect to (a/b) we refer to || |j] and for a discussion of 



the relationship of LLG to Landau gauge in perturbation theory we refer to [18]. To compute 



the lowest eigenmode of the covariant Laplacian we use the Arnoldi package ARPACK 



3 Comparison in the Coulomb phase 

We generate 1500 configurations on a 16 4 lattice at (3 = 1.05, fix them with the Laplacian and 
the modified iterative prescription and compare the results for the photon propagator. 

A gauge-invariant way to establish, for (3 > (3 C , that the photon is massless is to measure 
the (exponential) decay of the correlator between two planes of spatial Polyakov loops with a 
few nonzero momenta in the spatial directions. At large enough Euclidean time separation t 
one determines the effective mass m e g. After introducing the lattice momentum 

p fl = - sm(— ) = - sin(— ) (17) 
a 2 a 

and plotting m e s as a function of p 2 and fitting with 



m efl = ^ m 2 +p 2 , (18) 

the data ought to be consistent with m 2 = 0. Since the available (ordinary) momenta are 
p i = 27rnj/L i and good agreement with the dispersion relation (|18D typically requires p 2 ^Cl, 
the spatial extension of the lattice needs to be large. On our 16 4 lattice, results are acceptable 
for the two lowest momenta, |p| = l-27r/16 and |p| = \/2-27r/16, as displayed in Fig. 1 (top). We 
fit the correlator with a single and a double cosh and check whether the two effective masses 
are consistent with ([18]) where the mass has been set to zero, which indeed they are. 



Following Nakamura and Plewnia [11], we repeat this in a gauge-dependent way, using our 
two gauge-fixed ensembles. To that aim we determine the decay properties of the transverse 
photon-photon correlator (^), again for the same spatial momenta and with indices p = v which 
point orthogonal to both p and the propagation direction (associated with t). The result is 
displayed in Fig. 1 (bottom). The important point is that the link-link correlator yields the 
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Figure 1: Top: Correlator between planes ofPolyakov loops with once (+) or v2 times the minimum 
spatial momentum at /3 = 1.05. From these one might determine two points of the photon dispersion 
relation (^), but this gauge-invariant estimator is rather noisy. The same information may be gotten, 
after gauge-fixing, with considerably reduced error bars from the photon-photon correlator with the 
respective momenta in the modified iterative (bottom left) or Laplacian Landau gauge (bottom right). 



same information about the photon dispersion relation as the Polyakov loop correlator, but 
with much higher accuracy. This is easy to understand, since the main effect of the gauge 
transformations is to "unwind" the Polyakov loop, i.e. to remove multiples of 2n which only 
contribute to the noise. 

The main technical point of the present paper is that one may get more detailed information 
about the performance of either fixing procedure by considering the complete correlator 

cm = Y,Mx)My)e- ip(x - y) = Mp)M-p) (19) 

x,y 

for arbitrary four-momentum p. Since in the continuum the photon propagator is transverse, 
it makes sense to split it (on the lattice) into a transverse and a longitudinal piece, viz. 

CAP) = (0r ~ ^)D± + ^|| (20) 

with 

D^P 2 ) = l(tr(C ap )- P ^^) , D^) = P -^fR. (21) 

If the propagator is transverse (i.e. -D||(j3 2 ) = 0), the transverse coefficient D±(p 2 ) is simply 
obtained as 1/3 of the trace of the 
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Figure 2: Transverse component of the photon propagator at [3 = 1.05 on a 16 4 lattice, after modified 
iterative (left) and Laplacian (right) Landau gauge fixing. Since perturbatively D±(p 2 ) ~ 1/p 2 , the 
graphs p 2 D±(p 2 ) vs. p 2 (top) and 1/D±(p 2 ) vs. p 2 (bottom) are supposed to show an almost linear 
dependence. The modified ILG shows good rotational invariance, while the LLG shows a huge scatter 
of the data at a given p 2 (+). After the cylindrical momentum cut (x) the scatter on the r.h.s. is 
reduced (i.e. rotational invariance improves), but the IR enhancement in the Laplacian case persists. 



In Fig. 2 we display our results for the Coulomb phase. Since in the weak coupling limit 
D±(p 2 ) ~ 1/p 2 , we plot p 2 D±(p 2 ) vs. p 2 (top) and 1/D±(p 2 ) vs. p 2 (bottom) - both of them 
after modified iterative (left) and Laplacian (right) gauge-fixing. The striking fact is the huge 
amount of scatter in the data for the Laplacian case - data for a given p 2 do not collapse 
on a universal line, i.e. rotational invariance is seriously broken after Laplacian fixing in the 
Coulomb phase. A simple way to elaborate on this point is to apply a cylindrical momentum 
cut p0 |, i.e. to decompose each momentum according to p = p\\p +pi_ where p = ~(1, 1, 1, 1) 



and to retain only those momenta with a component orthogonal to (1, 1, 1, 1) that is smaller 
than some constant, e.g. |j3j_| 2 < 0.2. Those momenta out of the full set (+) which survive this 
cut are plotted with crosses (x). For large momenta the scatter is considerably reduced, but 
for the very low p 2 the cut gets less effective and the IR enhancement by which the Laplacian 
correlator differs from its modified iterative counterpart is unchanged. Towards the end of this 
article evidence will be presented that the Laplacian gauge suffers from DDS, and the present 
observation is consistent with that claim: If DDS persist in the Laplacian ensemble, then one 
would expect rotational symmetry to be broken (remember that DDS prefer to be planar) and 
expect low-lying momenta to be particularly affected, since DDS represent large-scale objects. 
In order to make contact with other studies, we fit the correlator to a functional form which 
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Figure 3: The longitudinal part of the correlator (fH|), multiplied by p 2 , in the Laplacian gauge at 
/3 = 1.05 (l.h.s.) and at /3 = 0.95 (r.h.s.). This confirms earlier findings that the LLG does not yield 
a purely transverse propagator, but the longitudinal admixture seems to be much more severe (in 
particular for small momenta) in the Coulomb phase (l.h.s.) than in the conhned phase (r.h.s.). 



allows for some IR enhancement; we use 

B <« 2 ) = Z q 2 il+a)+mW+al (22) 

with a "photon mass" m and an "anomalous dimension" a (apart from the renormalization 
factor Z). Then our results may be stated as follows: In the Coulomb phase the modified itera- 
tive prescription is consistent with a propagator ( f2^) with zero mass and vanishing anomalous 
dimension, whereas the Laplacian gauge suggests a positive anomalous dimension. 

An interesting point is to check the transversality of the propagator - we find it well 
obeyed in the modified ILG but seriously violated in the LLG. Fig. 3 shows ft 2 D\\ vs. p 2 for the 
Laplacian case both in the Coulomb phase (l.h.s.) and in the confined phase discussed below. 

To test whether the modified iterative and the Laplacian prescriptions differ on a deeper 
level (i.e. before the photon propagator is computed), we checked the distribution of link 
angles. They produce almost indistinguishible distributions (which are very close to a Gaussian 
centered about zero), i.e. they both comply with the perturbative preference for zero link angle. 



4 Comparison in the confined phase 

We generate 1500 configurations on a 16 4 lattice at /3 = 0.95, fix them with the Laplacian and 
the modified iterative prescription and compare the results for the correlator ([Tj|, |5|) . 

Following the order of investigation in the Coulomb phase, we start with the correlator 
between gauge-fields with 1 and \/2 times the minimum spatial momentum. The result is 
presented in Fig. 4. With either procedure the two correlators seem to lie on top of each other, 
indicating that the effective mass (|i"8|) is dominated by the m 2 contribution. This confirms our 
expectation that there is no massless excitation in the disordered phase, but the two graphs 
show different slopes so that the resulting effective mass is a gauge-dependent object. Moreover, 
it need not agree with what is usually called the photon mass, namely the asymptotic mass in 



the l + ~ plaquette-plaquette correlator IfLOl , which is a gauge- invariant quantity. The gauge- 
dependence of the effective mass of the correlator (^) for (3 < /3 C is also evident from Fig. 5 
and the Table below. Note finally that in the confined phase the Polyakov loop correlator is 
so noisy that in practice no signal can be extracted. 



8 




012345678 012345678 



Figure 4: Field (^) correlator with with once (+) or \/2 times the minimum spatial momentum at 
(3 = 0.95. This yields information about two points in the dispersion relation (|l^). In the confined 
phase the Polyakov loop correlator is so noisy that no signal can be extracted. The effective mass in 
the modified iterative gauge (l.h.s.) is smaller than that in the Laplacian gauge (r.h.s.). 



In order to investigate the difference between the two fixing procedures, we go through 
the same steps as before: We consider the complete correlator fllTf ) involving arbitrary four- 
momentum p^. The result for the confined phase is shown in Fig. 5. We plot p 2 D±(p 2 ) vs. p 2 
(top) and 1/D±(p 2 ) vs. p 2 (bottom) - both after modified iterative (left) and Laplacian (right) 
fixing. Unlike in the Coulomb phase, rotational invariance is not a problem, i.e. the data for a 
given p 2 collapse on a single line in either gauge. From a qualitative point of view it appears 
that in the confined phase fllTf ) is less sensitive to the details of the fixing procedure. On a 
quantitative level, however, the agreement is substantially less impressive. Again, we fit the 
correlators to the ansatz (|22|), which now seems much more adequate than in the ordered phase. 
The resulting parameters (cf. Table below) indicate that the confined phase is characterized by 
a distinctly nonzero mass parameter m and a positive anomalous dimension a, but the mass 
in LLG is substantially larger than that in the modified ILG. 

Again we check on the transversality. The modified iterative correlator is (almost perfectly) 
transverse, the Laplacian one is not. However, the non-transversality in LLG is seriously 
reduced compared to what we found in the Coulomb phase and rotational invariance is well 
respected in the transverse piece, too (see Fig. 3). In passing we note that these findings 
are consistent with the hypothesis (which will be discussed below) that the Laplacian Landau 
gauge lacks the ability to fully remove the DDS in the Coulomb phase - rotational invariance 
is restored in the confined phase, since the DDS percolate at the transition |17| . 

Finally, we check whether any difference between the two procedures may be seen in the 
distribution of the link angles, and here the answer is affirmative: In the confined phase, the 
distribution after Laplacian fixing deviates much more visibly from a Gaussian than that after 
iterative fixing. 



/3 = 0.95 

(confined) 


modif. ILG 

Z m a 


LLG 

Z m a 


all momenta 
cylind. cut 


5.25±0.04 1.10±0.02 0.41±0.01 
5.01±0.02 1.02±0.01 0.28±0.01 


4.29±0.05 1.33±0.02 0.40±0.02 
4.29±0.02 1.30±0.01 0.29±0.01 



Table 1: Coefficients of the fits of p 2 D±(p 2 ) vs. p 2 to (^) in the confined phase - see Fig. 5. The 
two lines represent the result without and with cylindrical momentum cut (see text). Note that the 
main uncertainty comes from the choice of the fitting method, not from statistics. 



9 




1 1 ' ' 1 1 ' 1 1 1 ' 1 1 ' ' 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 4 0.5 1 1.5 2 2.5 3 3.5 4 



Figure 5: Transverse component of the correlator ( |T^ ) at (3 = 0.95 on a 16 4 lattice, after modified iter- 
ative (left) and Laplacian (right) Landau gauge hxing. We plot p 2 D±(p 2 ) vs. p 2 (top) and 1/D±(p 2 ) 
vs. p 2 (bottom). Both gauges yield a propagator which is not perfectly rotationally invariant, since 
there is a certain amount of scatter among the data for a given p 2 (+). After the cylindrical mo- 
mentum cut (x) the scatter is much reduced. Note that the intercept, in the bottom line, of the 
correlator with the y-axis is larger in the Laplacian gauge than in the modified iterative case. 



5 Discussion and Conclusion 

Having established that on either side of the phase transition in the compact U(l) theory the 
correlator (|I9|) in the modified iterative Landau gauge (ILG) differs from that in the Laplacian 
Landau gauge (LLG), it is natural to ask which are the specific properties of the gauge fixing 
procedures which cause these differences. 

From the discussion of the gauge fixing procedures in sect. 2, it is clear that double Dirac 
sheets and zero-momentum modes represent the primary candidates for gauge artefacts which 
might spoil the properties of a gauge (cf. Hence, a natural thing to try is to omit the 

DDS suppression in the iterative procedure while keeping the step which reduces the ZMM 
background; this is what we have introduced as our "standard" ILG. The result of such a gauge 
fixing is presented in Figs. 6 and 7. In the Coulomb phase, the correlators with minimum or 
next-to-minimum spatial momentum (l.h.s. of Fig. 6) are flatter than those after modified ILG 
and LLG (cf. Fig. 1), indicating a negative m 2 . This means that the DDS tend to decrease the 
effective mass (|i~8|). In the confined phase the "standard" iterative correlators essentially agree 
with those in the modified iterative gauge. A more detailed picture is gotten by considering 
arbitrary four-momenta, as displayed in Fig. 7. The resulting longitudinal piece is essentially 
zero, hence only the transverse part is shown. The startling observation is that the "standard" 
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Figure 6: Field (^J correlator with with once (+) or y/2 times the minimum spatial momentum in 
the "standard" iterative Landau gauge at /3 = 1.05 (l.h.s.) and at /3 = 0.95 (r.h.s). The propagator in 
the Coulomb phase is flatter than its modified iterative or Laplacian counterparts at the same (3 value 
(cf. Fig. 1), while in the confined phase it is in good agreement with its modified ILG counterpart 
but not with the Laplacian one (cf. Fig. 4). 



ILG propagator shows the same amount of scatter and the same sort of IR enhancement in the 
Coulomb phase that we have seen in the case of the Laplacian correlator. On the other hand, 
in the confined phase the "standard" iterative propagator almost coincides with its modified 



counterpart. This is expected, since the DDS percolate at the phase transition [17]; hence the 
efforts of the modified ILG to remove them are futile. 

In an attempt to further substantiate the obvious suspicion that remnant DDS are respon- 
sible for the abnormal IR behaviour of the Laplacian propagator in the Coulomb phase and to 
shed light on the properties of the Laplacian gauge in the confined phase, we produce a scatter 
plot relating the value of the functional (0) to the remaining number of Dirac plaquettes in 
the gauge fixed configuration. The result is shown in Fig. 8. 

In the Coulomb phase the modified ILG always reaches the highest value of ([?]), among 
the three gauges considered. The "standard" ILG reaches (almost) the same extremum on 
the majority of configurations, but there is a considerable number of backgrounds on which 
it falls short quite drastically, and at the same time the total number of Dirac plaquettes is 
substantially increased. This correlation is paralleled in LLG, but there the amount by which 
it falls short (in either respect) has a continuous spectrum. In other words, the "standard" 
ILG differs from the LLG by the gap in the functional value which separates "success" from 
"failure". Hence it looks like only a limited number of configurations causes (simultaneously 
for the Laplacian and the "standard" iterative procedure) a problem both w.r.t. the gauge 
functional and the number of Dirac plaquettes. We checked that these are, in most cases, the 
same configurations for both gauges. On our lattice the DDS condition ([5p reads 2-16 2 = 512 
per plane. With a typical number of 0(2000) Dirac plaquettes at j3 = 1.05 in a modified 
ILG configuration, one would expect a configuration with a DDS sheet to have more than 
0(2500) Dirac plaquettes. This is precisely where the gap in the number of Dirac plaquettes 
for "standard" iterative fixing ends. Since the Laplacian fixing shows a continuous excess 
in the number of Dirac plaquettes, this correlation analysis cannot prove that DDS present 
a problem to the LLG procedure. However, the condition (||) being met, together with our 
previous observation that rotational invariance is seriously broken (in particular for small 
momenta, cf. Fig. 3) is strong circumstantial evidence. 

In the confined phase things look much simpler: Again, the modified ILG wins the com- 
petition for the extremization of the gauge functional (|7]). The "standard" ILG is in second 
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Figure 7: Photon correlator in the "standard" iterative prescription at (3 = 1.05 (left) and at (3 = 0.95 
(right). Each time the graphs for p 2 D±(p 2 ) vs. p 2 (top) and 1/D±(p 2 ) vs. p 2 (bottom) are shown. 
In the Coulomb phase the remaining DDS cause a massive IR enhancement which is absent in the 
modihed ILG counterpart but equally pronounced in LLG (cf. Fig. 2). In the confined phase the 
result agrees with the modihed ILG correlator (cf. Fig. 5). 

rank, but it is very close - both w.r.t. the functional value and the number of remnant Dirac 
plaquettes. The LLG is far off - again w.r.t. both criteria. The only surprise is by how much 
it misses the extremum and by how many Dirac plaquettes it deviates from the modified ILG. 
Again this establishes a link to our previous finding that the LLG respects rotational invari- 
ance in the confined phase: If unresolved DDS in the Coulomb phase are responsible for the 
shortcomings of the Laplacian gauge, then one would expect, since the DDS percolate at the 
transition (|T7|, that rotational invariance gets restored in the confined phase while the excess 
of Dirac plaquettes remains - and this is precisely what we see. 

While these preliminary observations clearly show a demand for future research, we would 
like to close with a few general remarks: 

• At first sight, it looks like some of our correlators are quite "fuzzy". Unlike many other 
groups we have refrained from averaging over momenta with the same p 2 ; if one would 
take the average D ± (p 2 ), then all propagators would become thin lines. However, as is 
obvious from the upper right plot in Fig. 2, such a procedure would only reveal a slight 
IR enhancement, and the information about the lack of rotational invariance would be 
lost. Hence we believe it is important to retain the full information about the spread of 
correlator points associated with a given p 2 . 
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Figure 8: Correlation between the final value of the functional ([?]) and the number of Dirac plaquettes 
in the fixed configuration in the Coulomb phase (l.h.s.) and in the conhned phase (r.h.s.). 



It is not clear whether the qualitative aspects of the propagators we have studied on a 
16 4 grid stay invariant under a change of the lattice size. On one hand the ability of the 
ILG to remove DDS through repeated attempts will degrade drastically on larger lattices; 
presumably the cut-off on the outer loop (as discussed in sect. 2) should be increased 
exponentially with the volume. On the other hand DDS on a larger lattice involve more 
Dirac plaquettes and become less frequent by purely entropic suppression. Which effect 
wins is unclear. Preliminary studies on a 32 4 grid indicate that the Laplacian correlator 
might be much better behaved on that lattice size - unfortunately the available statistics 
has prevented us from reaching any conclusion. Still, our 16 4 lattice is much larger than 



those which have been used in previous investigations of the modified ILG |14 



• After dealing with so many details of which quantity correlates with which other one, 
we wish to emphasize a simple point: We have presented evidence that all three gauges 
considered yield smooth configurations in either phase. At (3 = 1.05 we see no differences 
in the distribution of link angles at all, whereas at /3 = 0.95 the configurations tend to 
be slightly smoother with both types of iterative fixing than in the Laplacian gauge. 
Notably, in spite of the distribution of link angles being almost indistinguishable, the 
propagators ([19]) on these ensembles may be strikingly different. 

• The lack of several gauge fixing procedures to reproduce the known perturbative be- 
haviour in the Coulomb phase of the compact U(l) theory raises concerns about the 
validity of results for SU(2) or 577(3) in either standard iterative or Laplacian gauge. 
If already a Z7 (1) subgroup may cause a problem, one does not feel inclined to blindly 
trust the result of any fixing procedure in the non-Abelian case. The mass which we have 
extracted from the gauge-fixed field correlator is not gauge-invariant (see Table 1). This 
feature is likely to persist in the non-Abelian case below T c . Therefore, the quantitative 
meaning of gauge-fixed gluon mass measurements should be considered with caution. 
It will be interesting to see whether our observation that the deconfined phase is more 
sensitive to the details of the fixing procedure carries over to the non-Abelian case, too. 
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